Compact object coalescence rate estimation from short gamma-ray burst observations 
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Recent observational and theoretical work increase the confidence that short-duration gamma-ray bursts are 
created by the coalescence of compact objects, like neutron stars and/or black holes. From the observation of 
short-duration gamma-ray bursts with know distances it is possible to infer their rate in the local universe, and 
draw conclusions for the rate of compact binary coalescences. Although the sample of such events with reliable 
redshift measurements is very small, we try to model the distribution with a luminosity function and a rate func- 
tion. The analysis performed with a sample of 15 short gamma- ray bursts yields a range for the merger rate of 
75 to 660 Gpc^'^yr^^, with a median rate of 180 Gpc^^yr^^. This result is in general agreement with similar 
investigations using gamma-ray burst observations. Furthermore, we estimate the number of coincident obser- 
vations of gravitational wave signals with short gamma-ray bursts in the advanced detector era. Assuming each 
short gamma-ray burst is created by a double neutron star merger, the expected rate of coincident observations 
is 0.1 to 1.1 per year, when assuming each short gamma-ray burst is created by a merger of a neutron star and a 
black hole, this rate becomes 0.4 to 4.0 per year. 



I. INTRODUCTION 

Short Gamma-Ray Bursts (SGRB) are some of the most pow- 
erful explosions detected in the universe. They are intensive 
bursts of high energy gamma rays with a duration shorter than 
two seconds 1 1 1, with isotropic energies exceeding 10^^ ergs. 
The most common interpretation of their progenitors is a sys- 
tem of two compact objects, either two neutron stars (NS-NS) 
or a neutron star and a black hole (NS-BH), which undergo a 
merger IJllH. According to General Relativity the orbits of 
these objects shrink due to emission of gravitational waves, 
before they coalesce to form a single black hole. Because of 
the emission of gravitational waves, these objects are one of 
the primary sources for ground-based gravitational-wave de- 
tectors, such as LIGO |4| and Virgo |5 j. A detection of a GW 
signal on its own is already very interesting, as it allows to 
draw conclusions on the NS equation of state [6^ JJ, will al- 
low tests of general relativity in the strong-field regime jS), to 
set limits on the graviton mass ll9l [T0ll . and to test the Lorentz 
Invariance principle [llj. Detection of GW signals in coinci- 
dence with optical counterparts, such as a GRB, will not only 
increase the confidence that the GW signal is real, but also 
adds much more to the scientific outcomes. If, for example, 
the redshift z of the GRB can be measured, it becomes pos- 
sible to measure cosmological distances without the cosmo- 
logical distance ladder (see e.g. llT2l and references therein). 
Even if e.g. the Hubble parameter can be constrained from 
GW observations alone Itl3i . coincident detections will help 
to vastly improve the precision of such measurements lfT4l . 
GW-GRB observations will therefore become very important 
in the future in the area of cosmology. 

In this paper we estimate the local merger rate of compact ob- 
jects, as well as the number of expected GW-SGRB observa- 
tions in the advanced GW-detector era. Short GRB data with 
redshifts from SWIFT observations are fitted to a theoretical 
model. In order to avoid selective bias, only SWIFT GRBs are 



used, and each redshift association is carefully examined, re- 
jecting any GRB with an uncertain association. Furthermore, 
in contrast to a previous study |15 |, the individual measured 
luminosities of each GRB are incorporated to improve the re- 
sult, and a more complete data sample is used. 
This paper is organized as follows: Section [lljdesciibes the 
available data of SGRB used in this work and the fitting mod- 



els. Results are discussed in section III and compared to rate 
estimates from other works or obtained with alternative meth- 
ods. 



II. DATA AND MODELS 

We use the set of SGRB with reliable redshift measurements 
as listed in (15], complemented by three additional SWIFT 
SGRBs with reliable redshift measurements |16|. See Ap- 
pendix|A]for details on the selection of the additional GRBs. 
Restricting the analysis to 17 SWIFT GRBs with a reliable 
redshift measurement, we eliminate any instrumental bias 
with respect to GRBs detected by other GRB missions. 
The (apparent) luminosity of these GRBs is computed using 
the redshift and fluence information, as given in Table|l] To ap- 
proximate the flux F in the observer frame, which is the rele- 
vant frame for the detection threshold of the satellite, the GRB 
fluence S is divided by its duration. As the observed fluence 
depends on the spectral properties of the source and the energy 
response of the detector, 15 to 150 keV for the BAT instrument 
onboard SWIFT iflTl . and we observe sources over cosmolog- 
ical distances, we expect a spectral shift and a change of flu- 
ence for two identical sources at different distances. As the 
observed photon number spectrum can be well expressed by 
the Band function ifTSl . it is possible to calculate the spectral 
shift as a function of redshift (see, e.g., Refs |19 20|). How- 
ever, as most sources for the short GRBs are closer than z =1, 
it can be shown that, using the Band function with reasonable 
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FIG. 1. Distribution of the used SGRB as a function of luminos- 
ity and redsliift. Only 17 of these 19 SGRBs are used in this paper 
(unshaded area), with a redshift smaller than 1. 

parameters as given, e.g., in Ref. ||20ll , the effect of the spec- 
tral shift is less than 10%. As this error is small compared 
to other statistical and systematical errors, we will neglect the 
effect of spectral shift in this work. 
The apparent luminosity is 



L^Andl{z)F^4TTdl{z) — , (1) 

J 90 

where di,{z) is the luminosity distance for a given redshift 
z, S is the measured fluence and Tgo is the time over which 
the burst emits 90% of its total energy. We adopt a standard 
flat-A cosmology with Ho ~ 71 km s^^, JIm = 0.27 and 
r^A = 0.73. The approximate luminosity in the source frame 
is obtained by multiplying by (1 + z). 

Fig. [T]shows the distribution of the observed luminosities as a 
function of redshift, with the solid line indicating the approx- 
imate threshold on the detector's sensitivity: 

where we use the flux threshold as of Ref. 1*201 with i^iim = 
5 X 10^^ erg/s/cm^. Fig.[l]shows a general trend of missing 
GRBs above a redshift of z ^1; this is the so-called redshift 
desert [21]. It is mainly attributed to obscuration effects and 
red-shifted emission lines in the range between z ~ 1.1 and 
z ~ 2.2 ll22l . Restricting to GRBs with redshift smaller than 
1 leaves 15 data-points for this analysis. 
The data are fitted to a luminosity function, and several com- 
monly used luminosity functions were considered (e.g. a 
power law, the Schechter function etc., see Refs. fT5ll24ll25 l). 
Only the lognormal function produced a reliable fit to the data, 
so we use this function in the remainder of this work. This 
lognormal luminosity function is 




■ ■ Hopkins 
— Wilkins 

0.0 0.5 1,0 1.5 2.0 

redshift 

FIG. 2. Comparison of the used rate functions, which are explained 
in detail in Appendix [B] Although we consider only GRBs with red- 
shift 1, the rate function varies significant over this distance range. 
Not shown in this plot are the delay rate functions. 



cj){L) (X -exp I ^ 1 , (3) 

where is the mean (peak) value of the luminosity and a is 
the width of the distribution. These two parameters are deter- 
mined by fitting the function to the 15 GRBs in our sample. 
No evolutionary effects of the luminosity function are taken 
into account in this work, although some processes might de- 
pend on the metallicity of the progenitors 1.26. 27j. which is 
a function of redshift. Such effects may be included in future 
work. As this fit only utilizes the observed GRBs, we need to 
rescale the observed luminosity function 4)'{L) by the volume 
the satellite is sensitive to, assuming an isotropic distribution 
of SGRBs: 

m « , (4) 

where d^^^^{L) is the maximum luminosity distance to which 
a GRB can be seen (see Eq. Q). 

With the real luminosity function in place, it is now possible 
to calculate the number of observable SGRB within some red- 
shift distance z: 

Jo 1 + ^ dz' 7L,„i„(2') 

where dV{z')/dz' is the comoving volume element and A^o 
is a normalization factor The rate function R{z) describes 
the formation rate of binary systems per volume as a func- 
tion of redshift. Since a binary system of compact objects 
must be formed from massive progenitor stars, R{z) can be 
assumed to follow the star formation rate. Explicit expres- 
sions and references for R{z) used in this work can be found 
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TABLE I. List of SGRBs with a reliable redsiiift found by SWIFT between 2004 and 2011. The table includes measured flu- 
ences, and redshifts, and the luminosities derived from these quantities as described in the text. The redshift data are taken from 
1151 or listed in Appendix |A| while the fluence data are taken from 1 23 1 or taken from the Gamma-Ray burst Coordinated Network, 
http : //gen . gsf c ■ nasaTgov/gcn3-archive ■ html| ^GCN6623, ^GCN7148, ^GCN7761, "GCNSIS?, ^GCN9337, ®GCN10338, 
^GCN 11111,'*GCN11467. 



(jKB 


Duration 


Redshift 


Fluence 




LsJ 




[iU erg/cm J 


050416 


2.0 


0.6535 


3.7±0.4 


050724 


3.0 


0.2576 


10.0±1.2 


051016B 


4.0 


0.9364 


1.7±0.2 


051221 


1.4 


0.5465 


11.5±0.4 


060502B 


0.09 


0.287 


0.4±0.1 


060801 


0.5 


1.131 


0.8±0.1 


061006 


130 


0.4377 


14.2±1.4 


061201 


0.8 


0.111 


3.3±0.3 


070429B 


0.5 


0.9023 


0.6±0.1 


0707 14B 


64 


0.9225 


5.1±0.3' 


071227 


1.8 


0.384 


2.2±0.3^ 


080520 


2.8 


1.545 


0.6±0.1^ 


080905 


1.0 


0.1218 


1.4±0.2'' 


090510 


0.3 


0.903 


3.4±0.4^ 


100117 


0.3 


0.92 


0.9±0.1* 


100816 


2.9 


0.8049 


20.0±1.0'' 


101219 


0.6 


0.718 


4.6±0.3'* 



in Appendix [B] However, since the time between the forma- 
tion of the compact objects and the coalescence of the binary 
are presumably on the order of Gyr |28|, a significant delay 
with respect to the star formation rate can be assumed. This 
is taken into account by using a delayed rate function in Eq. |5] 
for two cases (see details in Appendix [B). Figure |2] shows the 
rate functions which differ significantly even over a range of 
redshifts smaller 1. This indicates, that it is very important to 
consider various theoretical rate functions in a work like this. 
Eq. (j5]l describes the observable number of SGRB, so the 
resulting number N'{z) need to be divided by a factor fn, 
describing the number of used SGRB to the total number of 
SGRB during the considered time span (2005 to 2011). Dur- 
ing that time span, 49 SGRBs were observed, of which 17 
constitute our sample. Therefore //? — 17/49, with an esti- 
mated error of about 20% (Vl7/49). The approximate time 
span spanning the range of data is about T = 6 yr We fur- 
ther neglect any downtime of SWIFT due to technical issues 
or other constraints, since this can be neglected compared 
to other sources of errors. Furthermore, SWIFT has a field 
of view of about 1.4 sr |17|, which means only a fraction 
Jfov —10% of the whole sky can be seen at a given time. 
Observations have shown that a short GRB could be cre- 
ated by a soft-gamma repeater, which makes up about 15% 
of all SGRB |I31|29l. We therefore assume that - 85% of 
short GRBs are created by the merger of two compact ob- 
jects, yielding a correction factor of /sgr = 0.85. SGRs are 
in general less luminous than ordinary SGRB, so events at 
larger redshifts might be composed purely out of SGRB. By 



using /sgr — 0.85, we will obtain a conservative limit on the 
merger rate. 

We also implicitly assume that each merger of two compact 
objects produces a GRB. In order to create a relativistic out- 
flow and therefore a GRB, a torus must be created around the 
newly formed black hole. Simulation and theoretical analyzes 
show that the formation of a GRB depends on several param- 
eters such as spin of the black hole, mass ratio or compact- 
ness of the neutron star (see, e.g., fT 'SO'l). Those parameters 
are hard to generalize, and assuming each merger produces a 
GRB is a conservative statement regarding the merger rates 
we are trying to estimate. 

SGRB and GRBs are believed to emit their radiation in a col- 
limated cone, whose half-opening angle 6 defines the fraction 
of the sky where the burst can be seen. This fraction is given 
by /b = 1 — COS0. The angle 9 is highly uncertain, espe- 
cially in the case of SGRB, as it depends on the model and the 
Lorentz factor of the outflow [3 1 1 . In the following, we adopt 
the value 1//^ = 15 ||32]| . corresponding to a half-opening 
angle 6* ~ 20° . 

With all these corrections in place, the total merger rate 
^merger {z) Can bc written as 



i?mcrgcr(2) = ^ . N {z) . (6) 

J Jb JFOV JR 

The expected rate of SGRB observations i?SGRB(^) is 
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Rsgrb{z) = ^Niz) . (7) 



III. RESULTS AND DISCUSSION 

The fit of eq. ^ to the data yield an observed peak luminosity 
of 4.2 X 10^2 ei-gs/s with a width of 2.26 x lO^^ grgs/s. The 
true peak luminosity can be found with Eq. Q yielding 1.3 x 
10*^ ergs/s. 
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The results of this analysis are summarized in Table III which 
show the approximated local coalescence rates of binary com- 
pact objects for the various rate functions. They are found 
to be between 150 and 330 Gpc~'^yr~^. With the large sys- 
tematic errors on the modeling and the assumptions made 
above, a more realistic range can be obtained if the minimum 
rate is halved and the maximum rate is doubled, yielding a 
rate between 75 - 660 Gpc~'^yr~^. The median rate is 180 
Gpc^'^yr^^. It should be noted that the upper limit on this 
rate comes from the model with unity rate function, not as- 
suming any evolution on the star formation over cosmological 
distances. The calculations have also been repeated with de- 
layed rate functions with delay times of 20 Myr and 100 Myr 
This resulted in a rate about 50% larger than without the de- 
lay, so the upper limit on our rate estimate need to be increased 
by 50%. Our final estimation for the merger rate in the local 
universe is thus 75 - 1000 Gpc^^yr^^. However, since 15 % 
of the observed SGRB could be SGRs instead, the real rates 
could be even larger. 

This result strongly depends on the half-opening angle of the 
SGRB jets. Figure [3] depicts the dependence of the rate on 
the opening angle 6. The smaller the angle, the larger the real 
number of merger events must be, since we have to be in the 
outflow cone to detect the SGRB in the first place. Assuming 
an half-opening angle of 10°, the merger rate could be as high 
as several thousand Mpc^-^yr^^. A more isotropic large half- 
opening angle of 60° yields a rate below 100 Mpc^'^yr^^. 
We also predict the number of coincident observations of 
SGRB with gravitational wave counterparts, assuming that a 
satellite with a field-of-view comparable to the field of view 
of SWIFT is operating at the time of advanced GW detec- 
tors. Using eq. (|7]i and the expected ranges for advanced 
LIGOA^irgo |33|, one expects 0.1 to 1.1 coincident observa- 
tions per year, if SGRB are created by a NS-NS merger (with 
a range of ^450 Mpc), and 0.4 to 4.0 coincident observations 



TABLE II. Merger rates summarized by neutron star observation and 
population synthesis computations, taken from ||33J . -Riow is a pes- 
simistic estimate, 7?re a realistic estimate and i?high an optimistic 
estimate. The results from this work are also included, ranging from 
75 to 1000 yr^^Gpc"^ 



source 


Rlow 


Rro 


Rhigh 


NS-NS (yr^^Gpc-^) 


10 


1000 


10000 


BN-NS (yr^^Gpc^^) 


0.6 


30 


1000 


this work (yr^^Gpc""^) 


75 (min) 


180 (median) 


1000 (max) 



10' ^ 



u 

o 



10' 

CO 

o 

^ 1 
10^ 
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FIG. 3. Coalescense rate of compact objects as a function of the half- 
opening angle 6. The red line indicates the median result while the 
grey area covers the possible range of rates including the variations 
from the different fit models and the systematic errors. The vertical, 
dashed line, indicates the opening angle used for the results in this 
paper. 



per year assuming a NS-BH progenitor (for which the range is 
^930 Mpc). These values include the systematic uncertainties 
underlying the estimates, as explained in the previous section. 
Since the advanced detectors will operate at least some years, 
these results suggest that it will be unlikely to not observe at 
least one coincident SGRB-GW event. 

If the number of operational GRB satellites changes, these 
numbers could change significantly. Assuming Fermi in 
operation during the time of advanced detectors as well 
as the planned SVOM mission fSTl and Lobster[35|, the 
coincident GRB-GW detection rate could be three times as 
high. Ensuring that at least one GRB mission is operational 
at the time of advanced detectors will be crucial to gain 
valuable information, like redshift, host galaxy, the galaxy 



TABLE III. Coalescence rates from modeling the observed 15 
SGRB, and the estimated detections per year for advanced 
LIGOMrgo with a reach of 450 Mpc for NS-NS systems and 930 
Mpc for NS-BH systems, respectively. 



Rate model 


rate 


NS-NS 


NS-BH 




Gpc-^yr"^ 




yr-^ 


Unity 


333 


171 


1123 


Wilkins 


171 


64 


577 


Hopkins 


155 


58 


524 


Fardal 


179 


74 


605 


Cole 


149 


52 


501 


Hemquist 


259 


125 


873 


Porciani 


156 


60 


492 


Porciani, 20 Myr 


252 


117 


849 


Porciani, 100 Myr 


259 


120 


874 



star formation rate etc, in addition to the parameters of the 
gravitational wave. Even if only one GRB mission will be 
operational, with an observing time of several years, there 
could be a good chance to get a few observations. 

The results obtained in this work can be compared to results 
obtained elsewhere. Two approaches are usually taken to esti- 
mate the rate of merger of compact objects. The first is to in- 
fer the rates from observed pulsar observations ll36l . The other 
approach is to employ population synthesis models f37 |. Both 
approaches inherit large statistical or systematical errors. A 
recent review article |33| summarizes those results, conclud- 
ing the real rate of merger events is somewhere between 10 
and 10000 Gpc^'^yr^^. Compared to that rate, the estimation 
in this work is more conservative. 

In this approach, the spatial distribution of observed short 
GRB with known redshift is fitted against a luminosity model 
and/or a rate function, from which the merger rate can be de- 
duced. Several investigations have been made in the past using 
this approach. Reference [38|, for example, finds a range of 
8-30Gpc ^^yr^^ for SGRBs, although this result is based on 
5 samples only. The analysis in ifTSl yields a much higher rate 
of about 7800 Gpc^^yr^^ with a sample of 7 (15) SGRB, but 
this work neglects the individual GRB luminosities. Finally, 
a recent study uses a different approach focusing on single 
GRB observations 1 19 1. The maximum distance at which each 
SGRB could be detectable by SWIFT is calculated, and con- 
verted into a final local rate. This study finds a merger rate of 
0.2 to 230 Gpc~'^yr~^ which overlaps with our rate estima- 
tion. 

Despite the various approximations used, our value of ex- 
pected detected mergers by LIGO is broadly consistent with 
that based on extrapolations from observed binary pulsar in 
our galaxy and population synthesis model [33 1- The method 
used in our work seems to be one of the most promising, 
especially with growing statistics of further observations of 
SGRBs. 
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|http : / / gen . gsf c . nasa . gov/ gcn3-archice . html[ 

GRB 080905 This short GRB has an optical afterglow and a 
host galaxy |39 |, making this redshift association reliable. 

GRB 090426 Although an optical afterglow is detected 
for this short GRB (GCN9255), spectral features indicate a 
long-GRB nature of this GRB (GCN9296). For that reason 
this GRB is not considered in our sample. 

GRB 100117 This short GRB has an optical afterglow 
and a host galaxy 140 il , making this redshift association 
reliable. 

GRB 100724 An optical afterglow has been observed 
(GCN 10971), but its spectral features do not allow a precise 
discrimination between a short or long GRB (GCN10976). 
For that reason this GRB is not considered in our sample. 

GRB 100816 This short GRB has an optical afterglow 
and a host galaxy (GCNl 1 123), making this redshift associa- 
tion reliable. 

GRB 101219 The XRT error circle for this GRB contains 
only one galaxy, from which a redshift obtained (GCN11518). 
This redshift association is therefore considered reliable. 



Appendix B: Rate functions 

This Appendix describes the rate functions used in Eq. |5]in its 
explicit forms. H(z) is defined as 

H{z) - Ho V(l + ^r^M + . (Bl) 

a. Porciani This function is actually the SF2 formulation 
given in Ref. BTI 
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Appendix A: Detailed redshift discussion 

This section lists all short GRBs discovered after 2009 with a 
redshift measurement. For each GRB, the redshift association 
is scrutinized. Only cases for which an optical afterglow or 
an unambiguous host galaxy was found, are considered in 
this work. The numbers in parenthesis correspond to specific 
circulars of the Gamma-Ray burst Coordinated Network, 



R{z) oc 



exp 3.4 z 
exp3.4z + 22 



b. Hernquist 

R{z) cx 



(B2) 



(B3) 



H{z) 



Ho 



2/3 



, a = 0.012 and /3 = 0.041 as taken 



with xiz) 
from Ref. 

c. Fardal 



R(z) (X -—-H(z) (B4) 

with pi = 0.075,^ = 3.7, P3 = 0.84 and a = 1/(1 + z) as 
taken from Ref. 
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d. Cole 



R{z) « 



(B5) 



with different parameters by different authors, summarized in 
the Table below. 



Ref 


a 


b 


c 


d 


Cole l44l 


0.0166 


0.1848 


1.9474 


2.6316 


Hopkins |45l 


0.0170 


0.13 


3.3 


5.3 


Wilkins (46| 


0.014 


0.11 


1.4 


2.2 



e. Delayed functions Because As the merger rate may not 
follow directly the star-formation rate, since 
Since the time between the formation of the compact objects 
and the coalescence of the binary due to the emission of grav- 
itational waves are presumably on the order of Gyr, a signif- 
icant delay with respect to the star formation rate can be as- 
sumed. Therefore, a delayed rate function is defined: 



Rt{z) 



At— R{Zret)P{t) (B6) 

i + ^/ 



with P{t) representing the probability distribution of the de- 
lay time. Population synthesis models Ii47i |48J suggest this 
distribution to be a power law 



P{t) cx t" 



(B7) 



with a ~ — 1 for t > itmin- z^ct is the retarded redshift, i.e. 
the redshift at the time the compact objects have formed. If 
T{z) is the lookback time at redshift z then 

^,et = (T(z) + t) , (B8) 
with the lookback time T{z) defined as (see, e.g., ||49|) 



V := fife(l + z)^ 
so the integral takes the form 



(CD 



3 /„(o) Vy/v + riA 

This has the generic solution 



(C2) 



T{z) = ^-^{Livo)~Liv,)) 



(C3) 



with setting 



(C4) 



This simple analytic expression can be used to calculate the 
lookback time for a given redshift, instead of evaluating the 
full integral. With this solution, it is now possible to derive an 
analytic inverse function of r(z). Solving Eq. (C3 i for L{yi) 
yields 



L{vi) = LM - 3^^^ = In E{T) 



(C5) 



The value of i(t;o) neither depends on z nor T, but is a mere 
function of fluj and JIa (see Eqs. (|C4[) and (|C1[)). Solving Eq. 



( C4 1 for z finally yields 



z{T) 



1/3 



i+E{T) y 

l-E{T)) 



1/3 



- 1 



(C6) 



T{z) = To 



dz' 



dz' 



(1 + + z'f + nk{i + z'Y + n^) 

(B9) 

The expression Eq. ( |B8| l contains an integration and function 
inversion, which is in general solvable by using numerical in- 
tegration methods. If constraining to flat cosmological models 
with Q,k = 0, it is possible to solve the integration and write 
down an an analytic expression for the lookback time. This 



helps to find an analytic expression for B8 and to reduce com- 
putational costs dramatically. The necessary steps are shown 
in Appendix [C] 



Appendix C: Analytic lookback times 



To solve Eq. B8 an inverse of the integral Eq. B9 is required. 
If restricting to flat cosmological models with = 0, one 
can substitute 
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